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A resonance model with magnetic connection for 3:2 HFQPO pairs 
in black hole binaries 
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ABSTRACT 

We apply epicyclic resonances to the magnetic connection (MC) of a black hole (BH) with 
a relativistic accretion disc, interpreting the high frequency quasi-periodic oscillations (HFQ- 
POs) with 3:2 pairs observed in three BH X-ray binaries. It turns out that the 3:2 HFQPO pairs 
are associated with the steep power-law states, and the severe damping can be overcome by 
transferring energy and angular momentum from a spinning BH to the inner disc in the MC 
process. 

Key words: accretion, accretion discs - black hole physics - magnetic fields - stars: individ- 
ual (GRO J1655-40, GRS 1915H-105, XTE J1550-564) - stars: oscillations - X-rays: stars 



1 INTRODUCTION 

As is well known, the high frequency Quasi-Periodic Oscillations 
, (HFQPO s) have been observed in several X-ray Binaries. As ar- 

■ gued by Ivan der Klij feOOd |2006|) , HFQPOs in X-ray binaries 
' probably originate from the inner edge of an accretion disc around 
, a black hole (BH) of stellar-mass, since millisecond is the natu- 
' ral timescale for accretion process in these regions. Although a 
I number of models have been proposed to explain the HFQPOs 

■ in X-ray binaries, no conse nsus has been reached on thei r phys- 
' ical origin (see a review by McCli ntock & Remillardl2006l , here- 
, after MR06). The 3:2 HFQPO pairs have been o bserved in a few 
. BH binaries, i.e., GRO J1655-40 (450, 300H z; iRemillard et al] 
• |l999l : IStrohmav ei"2Q0l': 'Remillard et al j|2002h, XTE J1550-564 
' (276, 184Hz; [Miller et al. 2001; Remill ard et al.ll2002h and GRS 

1915-^105 (168, 113Hz; MR06: IRemillard & McClintockl l2006l 
hereafter RM06). 

It has been pointed out that HFQPOs are generally asso- 
ciated with the steep power-law (SPL) state in BH X-ray bi- 
naries. Although the 3:2 HFQPO pa irs could be interpreted 
in some epicyclic resonance models i Abramowicz & Kluzniakl 
20011: lAbramowicz et alj l2003l : iKluzniak & Abramowiczl |2004 



Torok et al I2OO5I) , there remain serious uncertainties as to whether 



epicyclic resonance could overcome the severe damping forces and 
emit X-rays with sufficient amplitude and coherence to produce the 

HFQPOs (e.g., see a review in MR06). 

Very recently, iGan et al" I l l2009. hereafter G09l) proposed a 
model of magnetically induced disc-corona to interpret the spec- 
trum states in some BH binaries, in which the magnetic connection 
(MC) of a BH with its surrounding accretion disc can be produced 
naturally based on several minimal assumptions. It is found that the 



spectrum states in the BH binaries depend on accretion rate and BH 
spin, and henceforth this model is referred to as the MC model. 

In this paper we intend to combine the resonance model with 
the MC model to fit the association of the 3:2 HFQPO pairs with 
the SPL state in the BH binaries. This paper is organized as follows. 
In Section 2 we discuss the association of HFQPOs with the SPL 
state based on the MC model. In Section 3 we discuss the input of 
the electromagnetic energy onto the resonance modes to overcome 
the severe damping by transferring energy from the BH to the in- 
ner disc via the MC process. Finally, in section 4 we discuss some 
issues related to our model. Throughout this paper the geometric 
units G = c = 1 are used. 



2 ASSOCIATION OF HFQPOS WITH SPL STATES 

G09 resolved the dynamic equations of standard thin disc model, 
combining the MC effects in a typical disc-corona scenario, and 
the main points are summarized as follows. The interior viscous 
stress is assumed to be proportional to gas pressure, which is com- 
parable to magnetic pressure. The interior viscous process is domi- 
nantly governed by tangled small-scale magnetic fields. A fraction 
of the viscously dissipated energy is released in the disc, emitting 
eventually as blackbody radiation and supplying seed photons for 
Comptonization of corona, and the rest viscously dissipated energy 
is converted into energy in the corona, heating corona and main- 
taining its relativistic temperature. It is assumed that the large-scale 
magnetic field is generated from the small-sc ale magnetic field , be- 
tween which the relation is given as follows ( lLivio etalJll999h . 



Bo 



(h/r) ■ Bd. 



(1) 
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In equation (1) Bp is the poloidal component of the large-scale 
magnetic field anchored on the disc, and h is the half height of the 
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disc. The rotational energy of a BH can be extracted via the MC 
process, which contributes to the energy dissipation of the disc- 
corona system, and the mapping relation between the BH horizon 
and the inner disc can be determined based on the conservation of 
magnetic flux as follows. 



Be 




= 7r/2 



= ~Bh ■ 27rrH • 2M sin 6d6\r=r 



(2) 



where the left and right sides represent the magnetic fluxes thread- 
ing the disc and the BH horizon, respectively. The quantities M, 
rn and Bh are the BH mass, the horizon radius and the magnetic 
field at the horizon, respectively, and 9 is the angular coordinate of 
the horizon. The quantities A and D in equation (2) are relativistic 
correction factors in Kerr metric, and they read 



A = 1 + al/f'^ + 2a:/f=^ Z) = 1 - 2/f + al/f^ 



(3) 



where a, = J /M'^ is the BH spin defined in terms of the BH mass 
M and BH angular momentum J, and f = r /M is dimensionless 
disc radius. 

In our model we assume that the magnetic field is uniform 
at the horizon, and it is related to accretion rate Mo based on the 
equilibri um between mag netic pressure and the ram pressure given 



wee n mag i 
bv .Moderski et all ( 1 19971) . 

\/2MD/r2, 



(4) 



As argued in G09, the energy dissipation in the MC process 
can be derived by resolving relativistic conservation equations of 
energy and angular momentum for a disc with the MC. The X- 
ray spectra of the BH binaries can be simulated by resolving the 
magnetic disc-corona system with Monte-Carlo method^ 

According to the epicyclic resonance model jAschenbachl 
|2004 iTorok et al.ll2005h . the 3:2 HFQPO pairs are fitted by the 
radial and vertical resonance frequencies as follows. 



veji^r ~ 3/2 



3r ^ + 8a, f 



= -3/2 



where 



V4, 



o 2 — 2x1/2 

- 3a, r ) ' , 

2\l/2 



27rmBH(?''/^ + I 



(5) 
(6) 
(7) 



M^^ = 2.03 X lO^Hz. (8) 



In equations (5) — (8) is the Keplerian frequency, v,. and vg 
are respectively the radial and vertical resonance frequencies of the 
general relativistic disc, and rriBH is the BH mass in terms of the 
solar mass. 

By using equations (5) — (8) we obtain the resonance radii and 
the BH spins of the three BH binaries for the given BH mass in 
fitting the 3:2 HFQPO pairs. In addition, we can fit the SPL states 
of these sources by combining the resonance model with the MC 
model given in G09, and the concerned observational quantities and 
the fitting parameters are listed in Tables 1 and 2, respectively. 

The quantities ttibh, D and i in Table 1 are the BH mass in 
terms of the solar mass, the distance estimated of each source and 
the orbital inclination angle, respectively. The values of r32/rn,s in 
Table 2 are the resonance radii in terms of the radius of the in- 
nermost stable circular orbit (ISCO), being obtained by invoking 
vejvr — 3/2. The parameters m and a are the accretion rate in 
terms of Eddington accretion rate and the viscous efficiency, re- 
spectively. 
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Figure 1. Emerged spectrum of the three sources in the SPL state: (a) GRS 
1915+105 (1996 Nov. 28) fitted with mBH=18, »n=0.088 and a=0.39 in 
the upper panel, (b) XTE J1550-564 (1998 Sep. 18) fitted with »nBH=10.8, 
m=0.1 and «=0.391 in the middle panel, (c) GRO J1655-40 (1996 Aug.l) 
fitted with mBH=6.6, m=0.09, o=0.335 in the bottom panel. The total emis- 
sive spectrum and its thermal, comptonized and reflective components are 
plotted in thick-zigzag, thin-solid, dashed and thin-zigzag lines, respec- 
tively. The observation data (without error bars) of GRS 1915+105 are taken 
from Fig. 4.13 of MR06, and those of GRO J1655-40 and XTE J1550-564 
taken respectively from Figs. 1 and 2 of iRemillaid 12004.) . 
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Table 1. The observational quantities of tlie BH binaries witli the 3:2 HFQPO pairs. 



Source 




D (kpc)° 


i (degree) 


i'QPo(Hz)" 


i'X,SPL(i'Edd)'' 


GRS 1915+105 


10-18 


11 - 12 


70 ±2'= 


168,113 


1.1 


XTEJ1550-564 


8.4- 10.8 


5.3 ±2.3 




276,184 


0.5 


GRO J1655-40 


6.0-6.6 


3.2 ±0.2 


70"= 


450,300 


0.1 



""MROe; ' jpender et al.l <2004h : iGreineil j200lh : ferosz et al.l <2002l) ; ' Ifiailvn et alj J 19951) 



Table 2. The parameters in fitting the SPL state of the BH binaries with the 3:2 HFQPO pairs. 



Source 




a. 






m 


ct 


Photon Index 


-'^x(-f'Edd) 


-PdiscZ-Plotal 


i^PL/-Flotal 


GRS 1915+105 


10 


0.685 


1.931 


14.255 


0.250 


0.300 


2.875 


0.460 


0.274 


0.722 




18 


0.994 


2.895 


18.046 


0.088 


0.390 


2.466 


0.927 


0.383 


0.612 


XTEJ1550-564 


8.4 


0.888 


2.114 


19.670 


0.130 


0.300 


2.348 


0.442 


0.234 


0.757 




10.8 


0.990 


2.741 


16.913 


0.100 


0.391 


2.615 


0.897 


0.420 


0.576 


GRO J1655-40 


6.0 


0.955 


2.332 


19.463 


0.110 


0.300 


2.585 


0.549 


0.347 


0.646 




6.6 


0.989 


2.712 


19.380 


0.090 


0.335 


2.592 


0.754 


0.400 


0.593 



Inspecting Table 2, we find that the values of the concerned 
parameters are in accordance with the definition of the SPL state 
given in RM06, and some issues in the fittings are given as follows. 

(i) As shown in Table 2, the resonance radius r^i is less than 
the radius rom , i.e., rn is located in disc region covered by corona, 
and the radiation from the resonance mode is influenced directly by 
the corona. 

(ii) We take the upper limit to the BH mass with the lower 
limit to the distance from the observers to fit the spectra of the three 
sources, since the bigger BH mass and the less distance give rise to 
the stronger photon flux in the fittings based on the MC model. 

(iii) As shown in Table 2, the power-law component domi- 
nates significantly over the disc component for each source. The ad- 
justable parameters are only accretion rate and viscous coefficient 
a for the given 3:2 HFQPO pairs. It is found that the soft (<~10 
keV) and hard (>~10 keV) parts of the power-law component in 
the radiation increase and decrease with the increasing accretion 
rate, respectively, while the corresponding fractions increase and 
decrease with the decreasing viscous coefficient a. Thus, as shown 
in Table 2, we fit the bigger photon index in the SPL state by in- 
voking greater accretion rate. The high values of a given in Table 2 
arise probably from the differential rotation strengthened due to the 
transfer of angular momentum from the spinning BH to the inner 
disc in the MC process. 

(iv) The higher BH spins are required in fitting the 3:2 HFQPO 
pairs spin based on equations (5) — (8). On the other hand, accord- 
ing to the MC model, the higher spin corresponds to the lower hard- 
ness of the power-law component in the radiation, giving rise to 
the steeper photon index. Thus we have the association of the 3:2 
HFQPO pairs with the SPL state in these BH binaries. 

Based on the MC model we have the emerged spectra of SPL 
states of the three sources as shown in Figure 1 . 



3 ENERGY INPUT ONTO RESONANCE MODES 

In order to understand the input of the electromagnetic energy onto 
the resonance modes, we have an analogy of the BH magnetosphere 
to an electric circuit as shown in Figure 2. In both cases we find that 
the electromagnetic energy is transported to the load as Poynting 




(b) 

Figure 2. (a) A magnetosphere with closed field lines connecting a spinning 
BH with its surrounding disc is shown in the upper panel, where the thin and 
thick arrows represent the closed magnetic field and Poynting energy flux, 
respectively, (b) A steady circuit is shown in the lower panel, where the thin 
and dot-dashed arrows represent magnetic and electric fields, respectively. 
The thick solid and thin dotted arrows represent Poynting energy flux and 
the current in the circuit, respectively. 
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Table 3. Comparison of transporting electromagnetic energy in BH magnetosphere 
with DC circuit. 



Sources 



BH magnetosphere 



Steady DC circuit 



Extracted Energy 
Load 

Transport Direction 
Energy Conversion 
Electron Accelerating 
Radiation 
Energy Transport 



Spinning energy of a BH 
Inner disc 
From BH to inner disc 
Magnetic reconnection 
Induced Electric field 
Inverse-Compton 
Poynting flux 



Non-electromagnetic energy 
Resistance 
From battery to resistance 
Dissipation in Joule heating 
Steady Electric field 
None 
Poynting flux 



flux. As argued bv lWangetai] i l2003h the energy is extracted mag- 
netically from a fast-spinning BH, and transferred to the inner disc 
via the closed field lines in the MC process as shown in the upper 
panel of Figure 2. This scenario is very similar to the energy trans- 
fer from a battery to the resistive load in a circuit as shown in the 
lower panel of Figure 2. In the circuit electrons are accelerated by 
the electric field, resulting in the dissipation of the kinetic energy in 
Joule heating. In the MC process the electromagnetic energy trans- 
ferred into the inner disc could be converted into the kinetic energy 
of electrons by the induced electric field due to the magnetic recon- 
nection. In this way the severe damping in resonance modes could 
be avoided. 

On the other hand, the magnetic reconnection is more likely to 
occur due to the resonance of the plasmoids in the accretion flow, 
and the electrons in the corona could be accelerated to upscatter 
soft photons to hard photons in inverse-Compton radiation. This 
scenario is helpful to understand the association of the HFQPO 
with non-thermal radiation observed in the SPL state of the BH 
binaries. A rough comparison of the two processes of transporting 
electromagnetic energy is given in Table 3. 



4 DISCUSSION 

In this paper, we combine the epicyclic resonances in a relativistic 
accretion disc with the MC to interpret the 3:2 HFQPOs pairs in the 
BH X-ray binaries. It turns out that the 3:2 HFQPO pairs are asso- 
ciated with the SPL state, and the severe damping can be overcome 
by transferring energy from a spinning BH to the inner disc via the 
MC process. Some issues related to this scenario are discussed as 
follows. 



4.1 Weak Magnetic Field Assumption 

IlH j2002l) discussed the condition of quasi-steady state in the MC 
of a spinning BH with an accretion disc, and he assumed that the 
magnetic field is weak so that its influence on the dynamics of par- 
ticles in a thin Keplerian disc can be negligible. The "magnetic field 
assumption" is given by 

|V(BV87r)| (9) 

where p is the mass density of the disc and g is the gravitational 
acceleration produced by the BH. 

It is easy to check that the assumption given by equation (9) is 
valid in our scenario for the 3:2 HFQPO pairs, so that the epicyclic 
resonance radii can be determined in the context of general rela- 
tivity. Replacing the gravitational acceleration g in equation (9) by 
Q?r, we have 




10 20 30 40 50 

r/M 




Figure 3. The ratio of |V{i?^/87r)| to pQ?r versus the disc radius in 
the scenario for the 3:2 HFQPO pairs, (a) GRS 1915+105 with mBH=18, 
m=0.088 and q=0.39; (b) XTE J1550-564 with mBH=10.8, m=0.1 and 
0=0.391 and (c) GRO J1655-40 with mBH=6.6, m=0.09, a=0.335. 



RaUoJ^^^, (10) 

where Q. = 27ri/^ is the Keplerian angular velocity. Based on equa- 
tion (10) and the values of the parameters given in Table 1 we 
have the ratios versus disc radius as shown in Figure 3, from which 
we find that the ratios are much less than unity for these sources. 
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Thus the "weak magnetic field assumption" remains effective in 
our model. 

4.2 Estimating tlie BH Mass of H1743-322 based on 3:2 
HFQPO pair 

Since HFQPOs are probably produced in the inner disc region very 
close to ISCO, they might offer the most reliable measurement 
of BH spins (MR06; RM06). Compared with other methods for 
measuring spin, such as fitting the spectrum of the X-ray contin- 
uum, HFQPO method is much more simple, being independent of 
disc inclination relative to the BH's spin axis. Furthermore, the 3:2 
HFQPO pairs provide a very strict constraint to the BH spin, pro- 
vided the BH mass is estimated. 

A very narrow range of the BH spin, 0.989 < a* < 0.994, 
can be found in Table 2, which corresponds to the upper limit to the 
BH mass for fitting the 3:2 HFQPO pairs associated with the SPL 
states observed in the three BH binaries. Combining this narrow 
range of the BH spin with equations (5) — (7), we obtain a very 
simple relation between BH mass and the upper frequency of the 
3:2 HFQPO pairs as follows, 

mBHi'up = 3000, (11) 

where the upper frequency I'up = vq is equal to the vertical res- 
onance frequency given in equation (6). It is noted that mBH in 
equation (11) is the upper limit to the BH mass. Replacing v^^ in 
equation (11) by Sfo, we have almost the same relation /-'o(Hz) = 
931mBH given in MR06. 

The 3:2 HFQPO pair has been observed in the bright X-ray 
transient H1743-322. Although the mass of its BH primary has 
not been measured, its behavior resembles th e BH binaries XTE 
J1550-564 and GRO J1655-40 in many ways (Homan et alj |2005l : 
Kalemciet all [20 06: Remillard et al. 2002, 2006). Very recently, 
McClintock et al.l ((2009) pointed out that H1743-322 does contain 
a BH primary based on a detailed analysis of its strong similari- 
ties to XTE J1550-564. As a simple analysis, we have the upper 
limit to the BH mass of H1743-322, jtibh < 12.5, by substituting 
I'up = 240Hz into equation (11). We expect that this upper limit 
can be checked by future observations on H1743-322. 



ACKNOWLEDGMENTS 

This work is supported by the National Natural Science Founda- 
tion of China under grant 10873005, the Research Fund for the 
Doctoral Program of Higher Education under grant 200804870050 
and National Basic Research Program of China under grant 
2009CB 824800. We are very grateful to the anonymous referee for 
his (her) helpful suggestion for improving our work. 



Greiner 2001, Nature, 414, 522 
Homan J. et al., 2005, ApJ, 623, 383 

Kalemci E., Tomsick J. A., Rothschild R. E., Pottschmidt K., Cor- 
bel S., Kaaret P 2006, ApJ, 639, 340 
Kluzniak W., Abramowicz M. A., Lee W., 2004, AIPC, 714, 379 
Li L.-X., 2002, ApJ, 567, 463 

Livio M., Ogilvie G. I., Pringle J. E., 1999, ApJ, 512, 100 
McClintock J. E., Remillard R. A., 2006. In Compact Stellar X- 

ray Sources, ed. Lewin, van der Klis, pp. 157-214. Cambridge: 

Cambridge University Press. (MR06) 
McClintock J. E., Remillard R. A., Michael P R. et al. 2009, ApJ, 

698, 1398 
Miller J. M. et al. 2001, ApJ, 563, 928 

Miller J. M., Fabian A. C, Wijnands R., et al., 2002, ApJ, 570, 
L69 

Miller 2004, eprint: arXiv:astro-ph/0402I0I i'l 
Miller J. M., Fabian, A. C, Reynolds, C. S., et al., 2004, ApJ, 606, 
L131 

Miller J. M., Homan J., Steeghs D., et al., 2006, ApJ, 653, 525 
Moderski R., Sikora M., Lasota J. P. 1997, in Ostrowski M., 
Sikora M., Madejski G., Belgelman M., eds, Proc. International 
Conf., Relativistic Jets in AGNs. Krakow, p. 1 10 
Orosz J. A., Groot P J., van der Klis M., et al., 2002, ApJ, 568, 
845 

Remillard R. A., Morgan E. H., McClintock J. E., Bailyn C. D., 

Orosz J. A., 1999, ApJ, 522, 397 
Remillard R. A., Muno M. P, McClintock J. E., Orosz J. A., 2002, 

ApJ, 580, 1030 

Remillard R. A., 2004, in Kaaret P, Lamb F. K., Swank J. H., 
eds, Proc. AIP Conf. Vol. 714, X-Ray Timing 2003: Rossi and 
Beyond. Am. Inst. Phys., New York, p. 13 

Remillard R. A., McClintock J. E., 2006, ARA&A, 44, 49 (RM06) 

Remillard R. A., McClintock J. E., Orosz J. A., Levine A. M., 
2006, ApJ, 637, 1002 

Strohmayer T. E., 2001, ApJ, 552, L49 

Sobczak G. J., McClintock J. E., Remillard R. A., et al., 2000, 
ApJ, 544, 993 

Torok G., Abramowicz M. A., Kluzniak W., Stuchlik Z. 2005, 

A&A, 436, 1 
van der Klis M., 2000, ARA&A, 38, 717 

van der Klis M., 2006, in Lewin W. H. G., van der Klis M., eds. 
Compact Stellar X-Ray Sources. Cambridge Univ. Press, Cam- 
bridge, p. 39 

Wang D.-X., Ma R.-Y., Lei W.-H., Yao G.-Z., 2003, ApJ, 595, 109 



REFERENCES 

Abramowicz M. A., Kluzniak W., 2001, A&A, 374, L19 
Abramowicz M. A., Karas V., Kluzniak W., Lee W. H., Rebusco 

P, 2003, PASJ, 55,467 
Aschenbach 2004, A&A, 425, 1075 

Bailyn C. D., Orosz J. A., McClintock J. E., Remillard R. A., 

1995, Nature, 378, 157 
Fender R. P, Belloni T, Gallo E. 2004, MNRAS, 355, 1 105 
Gan Z.-M., Wang D.-X., Lei W.-H., 2009, MNRAS, 394, 2310 

(G09) 



